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Extending our previous work we classify the nonlinear magneto-optical response at low index 
surfaces of fee antiferromagnets, such as NiO. Antiferromagnetic bilayers are discussed here as models 
for the termination of bulk antiferromagnets. 
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Nonlinear optics has been proven to be very useful for 
the investigation of ferromagnetism at surfaces due to its 
enhanced sensitivity to twodimcnsional ferromagnetism 
. The magnetic effects are usually much stronger than 
in linear optics (rotations up to 90°, pronounced spin 
polarized quantum well state oscillations , magnetic 
contrasts close to 100%) l^-Q. Recently, Second Har- 
monic Generation (SHG) has been successfully applied 
to probe antiferromagnetism (visualization of bulk AF 
domains [|s|-|lO|]). The potential of SHG to study the sur- 
face antiferromagnetism has been announced in Rcf. ]Tl| ] 
extensively discussed in our previous paper |l2| . 

The practical importance of studies in this field fol- 
lows from the applications of antiferromagnetic (AF) ox- 
ide layers in devices such as those based on TMR (tun- 
neling magnetoresistance) , where a trilayer structure is 
commonly used. The central layer of TMR devices con- 
sists of an oxide sandwiched between a soft and a hard 
magnetic layer. For these technological applications it is 
necessary to develop a technique to study buried oxide 
interfaces. Such a technique can be SHG. One of the 
promising materials for the mentioned devices is NiO. 
However, to the best of our knowledge, the understand- 
ing of its detailed spin structure is scarce - even the spin 
orientation on the ferromagnctically ordered (111) sur- 
faces is not known. 

Our recent paper fl^l presented an extensive study of 
the nonlinear electrical susceptibility tensor xfi^ (* ne 
source for SHG within the electrical dipole approxima- 
tion), mostly for monolayer structures. It has been 
proven there that the spin structure of an antiferromag- 
netic monolayer can be detected by means of SHG. Also 
the possibility of antiferromagnetic surface domain imag- 
ing has been presented for the first time. As it was men- 
tioned in this previous work, bilayer spin structures are 
enough to account for the symmetry of a surface of a 
cubic antiferromagnet. Here we present an extention of 
that work to the (110) and (111) bilayer structures, thus 
completing our group theoretical analysis of low Miller- 
index antiferromagnetic surfaces. 



We follow exactly the group theoretical method de- 
scribed in Ref. At this point it is necessary to 
define the notions of "phase" and "configuration", used 
henceforth to classify our results. "Phase" describes the 
magnetic phase of the material, i.e. paramagnetic, ferro- 
magnetic, or AF. Secondly, the word "configuration" is 
reserved for the description of the magnetic ordering of 
the surface. It describes the various possibilities of the 
spin ordering, which are different in the sense of topol- 
ogy. We describe AF configurations, denoted by little 
letters a) to 1), as well as several ferromagnetic configu- 
rations, denoted as "ferrol" , "ferro2" , etc. The number 
of possible configurations varies depending on the surface 
orientation. All the analysis concerns collinear antiferro- 
magnets, with one easy axis. 

The tables show the allowed tensor elements for each 
configuration. The tables also contain the information on 
the parity of the nonvanishing tensor elements: the odd 
ones are printed in boldface. In some situations an even 
tensor element (shown in lightface) is equal to an odd 
element (shown in boldface), this means that this pair of 
tensor elements is equal in the domain which is depicted 
on the corresponding figure, but they are of opposite sign 
in the other domain. The parity of the elements has been 
checked in the operations 2 Z , 4 Z , and in the operation 
connecting mirror-domains to each other (for the defini- 
tion of the mirror-domain structure see Ref. |T^1 ). The 
domain operation(s) on which the parity depends is (are), 
if applicable, also displayed in the tables. If two or more 
domain operations have the same effect, we display all 
of them together. To make the tables shorter and more 
easily readable some domain operations (and the corre- 
sponding parity information for the tensor elements) arc 
not displayed, namely those that can be created by a su- 
perposition of the displayed domain operations. We also 
do not address the parity of tensor elements in the 6 Z nor 
3 Z operations for (111) surfaces nor any other operation 
that "splits" tensor elements, although these operations 
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also lead to a domain structure As has been dis- 

cussed in Ref. |Q it is possible to define a parity of the 
tensor elements for the 3 Z and 6 Z operations, however the 
tensor elements then undergo more complicated changes. 
The situations where the parity of the tensor elements 
is too complicated to be displayed in the tables are in- 
dicated by a hyphen in the column "domain operation" . 
For the paramagnetic phase, where no domains exist, we 
display the hyphen as well. For some configurations, none 
of the operations leads to a domain structure - in those 
configurations we display the information "one domain" . 
The reader is referred to Ref. p^] for the particularities 
of the parity check. 



A. (110) bilayer 

The previously described AF configurations of the 
(001) monolayer most commonly get split into two differ- 
ent configurations when a bilayer structure is considered. 
For the (110) bilayer it is not the case - only two of twelve 
AF configurations get split in this way, thus one obtains 
14 AF configurations of the (110) bilayer. Describing the 
results of our analysis we use the nomenclature of our 
previous article, i.e. the antiferromagnetic configurations 
are labeled by small letters. Only the four configurations 
that result from splitting of the two configurations of the 
monolayer structure are labeled by small letters with sub- 
scripts that carry the information about how they have 
been constructed from the (110) monolayer. For configu- 
rations with subscript "a" the lower layer is constructed 
by translation of the topmost layer by vector (0.5a, 0.5b), 
where a and b are interatomic distances within the (110) 
plane along x and y axes, respectively. For configura- 
tions with subsript "b" the vector of translation is (-0.5a, 
0.5b). This corresponds to the way we constructed the 
(001) bilayers in jjj. 

The configurations of the (110) monolayer structure 
are depicted in Fig. [I], and the way the bilayer is con- 
structed is depicted in Fig. ||. The tensor elements are 
presented in Table |. In general, we can observe five types 
of response. However, the possibility to distinguish AF 
configurations is not much improved compared to the 
(110) monolayer. Even the possibility to detect the mag- 
netic phase of the surface is not evident. 

As for the (001) surface jL2j, there is no difference in 
SHG signal between the monolayer and bilayer for the 
paramagnetic and ferromagnetic phases. For most AF 
configurations, however (confs. a), b), c), e), f a , ffc), g), 
h), j), k), and 1)) such a difference is present due to a 
lower symmetry of the bilayer. 



B. (Ill) bilayer 

In order to be consistent with our previous work [l^] 
we keep the same configuration names as in this earlier 
paper. That is why, for example, conf. b) is not present 
here. The spin configurations of the (111) bilayer are 
constructed from the configurations of the (111) surface 
of our previous work in the way that the spin structure 
in the second atomic layer is the same as in the topmost 
layer, but shifted accordingly to form a hep structure. 
Taking into account the spin structure of the second layer 
causes all the AF configurations to split, thus one obtains 
10 AF configurations of the (111) bilayer. The configura- 
tions are labeled by small letters (indicating their "par- 
ent" configuration) with subscript "a" if the mentioned 
shifting is along the positive x axis, and "b" if the shifting 
is along the negative S^y axis. 

The configurations of the (111) monolayer are depicted 
in Fig. |§| and the construction of the bilayer is depicted in 
Fig. |J. The corresponding tensor elements are displayed 
in Tab. O. The results are identical to those of our previ- 
ous work jL2j, where the second layer of the (111) surface 
was treated as nonmagnetic. This means that the spin 
spin structure of the second layer does not play any role 
for SHG, however the presence of the atoms in the second 
layer does. 

III. CONCLUSION 

From our results follows that SHG can probe maxi- 
mally two atomic layers of the surface of cubic two sub- 
lattice antiferromagnets, and only one of the paramag- 
netic or ferromagnetic surface. For the (111) surface, 
the spin structure of the second layer does not have any 
influence on SHG, i.e. it does not matter from the group- 
theoretical point of view if the investigated surface is a 
termination of a bulk antiferromagnet or a monolayer 
grown on a nonmagnetic substrate. However, these two 
situations can be very different from the band-theoretical 
point of view. 
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ther 3 Z nor 6 Z have this property - they must be applied 
at least 3 times to restore the original situation. Concern- 
ing the tensor elements, for 2 Z and mirror-operations a 
given tensor element may change its sign or not. Upon 
the application of the 4 Z operation it may become an- 
other tensor element (i.e. its x indices become y and vice 
versa) while its sign can change or remain unchanged. 
Upon the 3 Z operation, however, a single tensor element 
gets "split" into several tensor elements, and the notion 
of the conservation (or not) of the sign of the tensor ele- 
ment looses its sense. The tensor elements are described 
in the cartesian coordinate system, where quarters are 
the elementary entities. Rotation by an angle other than 
90° and its multiples cannot be described as interchang- 
ing the axes and a possible modification of their signs. 
However, it is possible to treat the parity of tensor ele- 
ments in the 3 Z and 6 Z operation in a more complicated 
way. 
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FIG. 1. Spin configurations of an fee (110) monolayer. Except for confs. "ferro3", g), h), and i), the arrows always indicate 
in-plane directions of the spins. In confs. "ferro3", g), h), and i) (<gi) denote spins pointing along the positive (negative) 
z-direction, respectively. 
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FIG. 2. Spin structure of an antiferromagnetic (110) bilaycr constructed from a shift of the monolayer, where two different 
shiftings are applied. Filled (empty) circles represent the topmost (second) layer. The rightmost panel shows the conventional 
unit cell for the resulting bilayer structure. Here, conf. a) of the (110) surface serves as an example. 
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FIG. 3. Spin configurations of an fee (111) monolayer. Except for confs. "ferro5", k), 1), and m), the arrows always indicate 
in-plane directions of the spins. In confs. "ferro5" and k) (<g>) denote spins pointing along the positive (negative) z-direction, 
respectively. 
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FIG. 4. Spin structure of an antiferromagnetic (111) bilayer constructed from a shift of the monolayer, where two different 
shiftings are applied. Filled (empty) circles represent the topmost (second) layer. Here, conf. f) of the (111) monolayer serves 
as an example. The rightmost panel displays the conventional unit cell for the resulting bilayer structure of conf. f a ). 
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TABLE I. Nonvanishing elements of xi? f° r an s P m configurations of the (110) surface of a fee lattice. We denote xYj 
by ijk. Odd elements are in bold if a parity operation exists. The configurations are depicted in Fig, [jj, 



(2") 



conf. 



point group symmetry operations domain operation non-vanishing tensor elements 



para 
ferrol 

ferro2 



mm2 
m 



ferro3 
ferro4 

AF: 

a) , g), j) 

b) ,h),k) 

c) , d), 1) 
e), fa), ft) 



1,2* 
1 



1,2* 
1,2, 

1,2, 
1,2* 
1,2, 
1 



2z , 2j;- 



2a:, 2 y 

2 Z 



io), id) 



mm2 



one domain 



xxz = xzx, yyz = yzy, zxx, zyy, zzz 
xzx = xxz, xxy = xyx, yxx, yyy, yzz, 
yyz = yzy, zxx, zyy, zzz, zyz = zzy 
All the elements are allowed: 

xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
xyz = xzy, xxz = xzx, yyz = yzy, yzx = yxz, 
zxx, zyy, zzz, zxy = zyx 
xxx, xyy, xzz, xxz = xzx, yyz = yzy, 
yyx = yxy, zxx, zyy, zzz, zzx = zxz 

xxx, xyy, xzz, xxz = xzx, yyz = yzy, 

yyx = yxy, zxx, zyy, zzz, zzx = zxz 

xzx = xxz, xxy = xyx, yxx, yyy, yzz, 

yyz = yzy, zxx, zyy, zzz, zyz = zzy 

xyz = xzy, xxz = xzx, yyz = yzy, yzx = yxz, 

zxx, zyy, zzz, zxy = zyx 

All the elements are allowed: 

xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx — zxz, zxy = zyx 
xxz = xzx, yyz = yzy, zxx, zyy, zzz 
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TABLE II. Nonvanishing elements of xL f° r an s P m configurations of the (111) surface of a fee lattice. More monolayers 
are taken into account. We denote Xijk by ^ ne configurations are depicted in Fig. ^. 

conf. point grp. symmetry ops. domain op. non-vanishing tensor elements 



para 
ferrol 



3m 1, ±3*, 2 y , 2s( X y), %s(,-xy) 



ferro3 
ferro5 



1: 2j/ 

1,±3 Z 



AF: 

&a), a&), i a ), 

it), ka), kb) 

Co), C 6 ), 
fa), f„) 



zxx = zyy, xxz = xzx = yyz = yzy, zzz, 
xxx = -xyy = -yxy = -yyx 
All the elements are allowed: 

xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
xxx, xyy, xzz, xxz = xzx, yyz = yzy, 
yyx = yxy, zxx, zyy, zzz, zzx = zxz 

xxx = -xyy = -yxy = -yyx, xyz = xzy = -yxz = -yzx, 
xzx = xxz = yyz = yzy, xxy = xyx = yxx = -yyy, 

zxx = zyy, zzz 



xxx, xyy, xzz, xxz = xzx, yyz = yzy, 
yyx = yxy, zxx, zyy, zzz, zzx = zxz 
All the elements are allowed: 

xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
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